Nonequilibrium energetics of a single Fi-ATPase molecule 
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Molecular motors drive mechanical motions utilizing the free energy liberated from chemical 
reactions such as ATP hydrolysis. Although it is essential to know the efficiency of this free energy 
transduction, it has been a challenge due to the system's microscopic scale. Here, we evaluate the 
single-molecule energetics of a rotary molecular motor, Fi-ATPase, by applying a recently derived 
nonequilibrium equality together with an electrorotation method. We show that the sum of the heat 
flow through the probe's rotational degree of freedom and the work against external load is almost 
equal to the free energy change per a single ATP hydrolysis under various conditions. This implies 
that Fi-ATPase works at an efficiency of nearly 100% in a thermally fluctuating environment. 

PACS numbers: 05.70.Ln,05.40.Jc,87.16.Nn 



Fi-ATPase, a water soluble part of F Q Fi-ATP syn- 
theses, is a rotary molecular motor [D-Q. The central 7 
shaft rotates unidirectionally within a cylinder consisting 
of three a and three ft subunits while hydrolyzing ATP to 
ADP and phosphate (Fie. [IK) . A single ATP hydrolysis 
triggers a 120° rotation [4j,|5|. By fixing an Fi motor to a 
glass surface and attaching a probe filament or particle to 
the 7 shaft, we can observe its ATP-driven rotations us- 
ing conventional optical microscopes at a single molecule 
level 0, 13) El- In cells, combined with the membrane em- 
bedded proton driven motor F Q , they couple ATP synthe- 
sis/hydrolysis and proton flow. Thus, Fi-ATPase plays 
a central role in biological energy transduction. There- 
fore, it is crucial to reveal its energetics, especially the 
efficiency, to understand the principle of its operation^]. 
Although Fi-ATPase is known to be highly efficient [4|, 
well-controlled precise evaluation has not been achieved. 
In this study, we evaluated the energetics of the Fi- 
ATPase by measuring thermodynamic quantities of its 
probe particle using a new nonequilibrium equality [§- 
llfj ] together with an electrorotatoin method developed 
for torque manipulation of microscopic ob iects |lll4l3| . 

In single molecule experiments, what we can access is 
only the probe attached to molecular motors since motor 
proteins are quite small with a dimension of around 10 
nm. Accordingly, a methodology to extract energetic in- 
formation of motors from the probe motion is required. 
Although the probe is much larger than the motor pro- 
tein, it is still sufficiently small that thermal fluctuations 
have a dominant effect on its behavior. The probe ex- 
changes energy with an environment through thermal 
fluctuations as heat 14|, [l5[ . Consider a microscopic ob- 



ject moving in a viscous fluid with a frictional coefficient 
T and a velocity v(t). It feels a force of —Tv(t) + 
from the fluid, where is the thermal fluctuating force. 
Then, the heat flow per unit time from the system to the 
heat bath is naturally defined as J = ([Tv(t) — v(t)), 
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FIG. 1: (Color online) a, Schematic of Fi-ATPase molecule. 
The central 7 shaft rotates unidirectionally when ATP is hy- 
drolyzed to ADP and phosphate. 5 and e subunits are omitted 
for simplicity, b, Typical trajectory at low ATP concentra- 
tion (0.4 fM ATP, 0.4 (M ADP, 1 mM Pi ). An ATP binding 
reaction triggers a 120° step. 



where (•) is the ensemble average 3, 13] ■ However, since 
is usually not accessible in experiments, measure- 
ment of J in small fluctuating systems has been diffi- 
cult. Recently, a nonequilibrium equality that enables us 
to calculate J from experimentally accessible quantities 
was derived by Harada and Sasa for Langevin systems 
in nonequilibrium steady states fil-Hcij] and experimentally 
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where v s is the mean velocity. C(f) is the fluctua- 
tion of the velocity around v s calculated as the Fourier 
transform of the self time correlation function C(r) = 
([v(t + t) — v s ][v(t) — v B ]). R(f) is the Fourier trans- 
form of the velocity linear-response function ; (f(t)) N = 



0(N 2 ), where (v(t)) N is the 



v * + f-JsR(t-s)N(s) 
ensemble average of v(t) under a sufficiently small probe 
torque N(t). R(f) indicates the sensitivity of the veloc- 
ity to a perturbation at /. R'(f) is the real part of R(f). 
T is the temperature. Boltzmann constant was set to one 
for the briefness. Around an equilibrium state, fluctua- 
tion response relation (FRR) C(f) = 2TR'(f) is known 
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to hold [2(|, meaning that a system with large fluctua- 
tions is highly sensitive to an external perturbation and 
vice versa. On the other hand. FRR is generally violated 
in nonequilibrium systems [21- 23]. Therefore, the magni- 
tude of the integral in (pj indicates how far the system 
is from equilibrium and vanishes around the equilibrium 
state. Thus, the first and second terms of the equality 
correspond to the dissipation originating in the steady 
motion and that originating in the nonequilibrium fluc- 
tuations, respectively. It should be noted that heat can 
dissipate through multiple degrees of freedom, including 
rotations, swings, radial fluctuations, etc., in a rotational 
system such as Fi-ATPase. The equality enables us to 
track the energy flow in terms of the amount of heat dis- 
sipation through a degree of freedom of our interest by 
measuring the fluctuation and response of that degree of 
freedom. In the following, we focus on the rotational de- 
gree of freedom of the probe, where v(t) is the rotational 
velocity and T the rotational frictional coefficient. 

The experimental setup including the torque calibra- 
tion was essentially the same as that in [l3| (see [3] for 
the details). Fi molecules derived from a thermophilic 
Bacillus PS3 with mutations (His 6 -aC193S/W463F, 
Hisio-A 7S107C/I210C) [| were fixed on a cover slip 
functionalized by Ni 2+ -NTA (Fig. EJa)). Rotations of 
a dimeric probe particle with a diameter of 0.287 /nm 
attached to the 7 shaft were observed at 24±1 °C in a 
buffer containing 5 mM MOPS-K, indicated amount of 
ATP, ADP, Pi, and 1 mM excess of MgCl 2 over ATP and 
ADP (pH 6.9) on a phase-contrast microscope (Olympus) 
equipped with a high-speed camera at 1800 Hz. Torque 
necessary for an external load and the response measure- 
ment was induced using an electrorotation methodflll- 
Fui ] as illustrated in Fig. [2j The magnitude of torque is 
constant in the meaning that it does not depend on the 
particle motion but only on the delay, and is proportional 
to the square of the amplitude of the applied voltage and 
the volume of the dielectric objects. To apply a sinusoidal 
torque for the response measurement, we modulated the 
amplitude of the applied 10 MHz AC voltage periodi- 
cally at various frequencies in the range from 0.3 Hz to 
600 Hz. Temperature rise due to the electric field was 
at most 0.3°C, which was smaller than the temperature 
variance of the environment. 

At a low ATP concentration, the probe attached to 
the 7 shaft of Fi-ATPase rotated unidirectionally with 
discrete 120° steps (Fig. [TJd). Intervals between steps 
correspond to the waiting time for ATP binding. On 
the other hand, some particles showed bidirectional rota- 
tional Brownian motion. They were supposed to be at- 
tached to collapsed Fi molecules or glass surface through 
some non-specific bindings. As a control, we measured 
C(f) and R'(f) of such fluctuating particles in a buffer 
lacking ATP, ADP, and P ; . Figure EJa) shows that the 
FRR holds to a good extent, or C(f) and 2TR'(f) coin- 
cide well, except in the high frequency region, where sys- 
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FIG. 2: (Color online) a, Experimental setup. A dimeric 
polystyre particle is attached to the 7 shaft of Fi-ATPase 
fixed on the top glass surface. Four electrodes are coated on 
the bottom glass surface. Distance between electrodes is 50 
fim . b, Electrorotation method. A rotating electric field at 
a frequency of 10 MHz was generated at the center of four 
electrodes by applying sinusoidal voltages with a phase shift 
of 7r/2. Dielectric objects in this rotating electric field have 
a dielectric moment rotating at 10 MHz. The phase delay 
between the electric field and the dielectric moment resulted 
in the torque on the objects. By periodically modulating the 
amplitude of the applied voltage, we obtained a sinusoidal 
torque. Not to scale. 



tern noise, aliasing effects due to the finite sampling rate, 
and the finite exposure time distort the data. On the 
other hand, for continuously rotating particles driven by 
Fi molecules in a buffer containing ATP, ADP, and Pi , 
the FRR was violated in the low frequency region (Fig. 
C2h>)). We suppose that the Lorentzian-like spectrum in 
the high frequency region reflects Brownian motion in 
ATP waiting dwells, and the low frequency region reflects 
stepwise rotation. Then, we examined the dependence of 
C(f) and R'(f) on ATP, ADP, and P ; concentrations. 
As ATP and ADP concentrations were simultaneously 
increased, average rotational velocity and both C(f) and 
R'(f) increased (Fig. [3b-d). At low ATP concentra- 
tions, since Fi-ATPase spends most of its time trapped in 
the ATP waiting state, the fluctuation and response are 
suppressed. On the other hand at high ATP concentra- 
tions, it behaves like a Brownian particle running down 
a potential gradient, causing the increase in C(f) and 
R'(f)- When we held the ATP concentration constant 
and decreased the ADP and Pi concentrations, the aver- 
age velocity was nearly the same. In this case, although 
R'(f) remained almost constant (Fig. [3fc,e), C(f) in- 
creased significantly. This is probably due to the change 
of fluctuations triggered by the binding and release of 
ADP and Pi , while an apparent difference in the rota- 
tional behavior was not evident. At last, we applied a 
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FIG. 3: (Color online) Typical examples of fluctuations C(f) 
(lines) and responses 2TR'(f) (circles) of the rotational veloc- 
ities. The shaded area corresponds to the half of the integral 
in The rotational trajectory is shown in the inset, rots : 
rotations, a, Rotational Brownian particle in the absence of 
ATP, ADP, and P; . b-f, Particles driven by Fi-ATPase at 
indicated concentrations of substrates, f, Under a constant 
opposing torque of 17.7 pN-nm/rad. 



constant external torque in the opposing direction of the 
ATP-driven rotations. A torque of 17.7 pN-nm/rad de- 
creased the mean rotational velocity of a probe from 5.08 
Hz to 2.58 Hz (Fig. [5f). Under such an opposing torque, 
the extent of FRR violation, or the discrepancy between 
C(f) and 2TR'(f), was less than that in the absence of 
the constant opposing torque (Fig. [3h, f). 

From the above results, the amount of heat dissipa- 
tion through the rotational degree of freedom was esti- 
mated using ([1]). The integral in ([1]) corresponds to the 
twice of the area between C(f) and 2TR'(f) shaded in 
Fig. |3] since C(f) and R'{f) are even functions. We 
limited the integration in the range up to 300 Hz to 
avoid the systematic noise in the high frequency region. 
The frictional coefficient T was calculated from the av- 
erage of C(f) around 300 Hz by assuming the FRR in 
high frequency regions as C(f) = 2T/T. The amount 
of heat dissipation per 120° rotation was calculated as 
<3rot = J/iv s = Q s + Qv, where Q s = Tv^/3v s and 
Q v = Tjdf \c(f) - 2TR'(f)} /3v s . Q s and Q v are sep- 
arately plotted in Fig. |4] ATP hydrolysis and mechan- 
ical motion are known to be tightly coupled in the hy- 
drolytic reactions of Fi-ATPase; a single ATP hydrolytic 
cycle corresponds to a single 120° stepji], Thus, the 



free energy consumption during 120° rotation is A/i = 
A^atp - A^adp - MP, - AAi + T In ( [ATP] / [ADP] [Pi] ) . In 
Fig. BJUiii), the ratio [ATP]/ [ADP] [P;] and accordingly 
A/z are held to be constant except for a minor difference 
due to changes in the ionic strength, concentration of 
free Mg 2+ , etc. Under these conditions, Q rot was almost 
the same and equivalent to lS.fi shown as a horizontal 
bar, although Q s and Q v varied depending on the ATP 
concentration. The increase of ATP concentration re- 
sulted in the increase of Q s and decrease of Q v . This 
is consistent with that the 7 shaft rotates rapidly and 
smoothly at high ATP concentrations. When we held 
the ATP concentration constant at different A/x, the av- 
erage velocity and accordingly Q s were almost the same 
(Fig. IDJii, v)). On the other hand, Q v increased to- 
gether with A/i. As a result, Q ro t was again comparable 
with A/j,. Under a constant external load N, Fi per- 
forms work W = N x 120° against the load during 120° 
rotation. Although Q s and Q v decreased simultaneously 
under an external torque (Fig. HJiv)), the sum W + Q TOt 
was almost equal to A/i. Thus, we obtained an energy 
balance relation W + Qrot — A/x under all the conditions 
we examined. In addition, the Stokes efficiency defined 
as the ratio (X/ (A/i — W) was always less than unity as 
expected [3, [23 (Fig. |4j. High Stokes efficiency implies 
that the driving force of a motor is nearly constant. 

Since we attached a large probe particle to the 7 shaft 
with a relatively soft elastic linker, the time scale of the 
7-shaft motion is much faster than that of the particle. 
As the 7 shaft rotates, it performs a work on the linker, 
which is stored in the linker as an elastic energy. This 
elastic energy then dissipates through the delayed motion 
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FIG. 4: (Color online) Amount of heat dissipation and work 
per 120° rotation. Gray-scale bars correspond to different 
molecules. Colored bars correspond to the means. The val- 
ues of the standard free energy difference A/i° are different 
among literatures [261-28] . The horizontal gray bars indicate 
the range of A/i calculated on the basis of values reported in 
literatures [24, 126H28I ] by a method developed by Krab and 
van Wezel [2g ; 72.4-78.4 pN-nm (i-iv) and 109-115 pN-nm 
(v). Error bars indicate the standard deviations. 
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of the particle as Qrot or used to oppose external load 
if subjected to external load. The work on the linker 
during 120° rotation can not exceed A/i since it is the 
work against the conservative elastic force. Therefore, 
W + Q rot is limited by A/i. Note that W + Qrot is an 
extractable work by coupling an appropriate system to 
the probe. Thus, our result W + Q ro t ~ A/i implies that 
Fi motor works at an efficiency of nearly 100%. When 
the time scale of the probe is faster than that of the 7 
shaft by attaching a sufficiently small probes with a stiff 
linker, it is expected that the particle follows the motion 
of the 7 shaft without delay. If this is the case, in addition 
to A/z, there may be inextractable heat absorption or 
dissipation through the rigid complex of the probe and 
shaft. In our system with a large probe, it is possible that 
there is such an additional heat flow through the rotation 
of the 7 shaft. However, this heat does not conduct to 
the particle and is directly exchanged between the shaft 
and surrounding medium. Thus, by insulating the probe 
from the 7 shaft, we measured only the extractable work 
of Fi-ATPase and found that it works at an efficiency of 
nearly 100%. 

In summary, we measured the violation of the fluctua- 
tion response relation of Fi-ATPase at a single molecule 
level using an electrorotation method. From the extent 
of this violation, we calculated the thermodynamic quan- 
tities of the probe using a nonequilibrium equality and 
found that Fi-ATPase works at an efficiency of nearly 
100% under a large variety of conditions. The free en- 
ergy conversion at 100% efficiency does not contradict 
the thermodynamic laws, yet it is surprising that such a 
highly efficient machinery exists and is actually working 
in cells. It is also worth noting that Fi-ATPase adapts 
to a large variety of conditoins to retain this efficiency. 
Since the rotation of the 7 shaft is the sole handle to 
transmit energy between Fi motor and F Q motor in cells, 
this high efficiency might have been developed for the ef- 
ficient coupling of ATP hydrolysis / synthesis and proton 
transport during a long evolutionary process. Including 
this possibility, detailed mechanism and meaning of such 
highly efficient energy transduction remain a subject for 
future study. 

We appreciate critical discussions with Masaki Sano, 
Takahiro Harada, Makito Miyazaki, Felix Ritort, Ken 
Sekimoto, and Kazuhiko Kinosita, Jr. We thank Shigeru 
Sugiyama for the technical advice on the electrorotation 
method. This work was supported by Japan Science and 
Technology Agency (JST) and Grant-in-Aid for Scien- 
tific Research on Priority Areas, 18074001, 17049015, 
19037022, 18031033 (to E.M.), and 21740291 (to S.T.). 



[1] P. D. Boyer, Biochim. Biophys. Acta 1140, 215 (1993). 
[2] J. P. Abrahams, A. G. W. Leslie, R. Lutter, and J. E. 

Walker, Nature 370, 621 (1994). 
[3] H. Noji, R. Yasuda, M. Yoshida, and K. Kinosita, Jr., 

Nature 386, 299 (1997). 
[4] R. Yasuda, H. Noji, K. Kinosita, Jr., and M. Yoshida, 

Cell 93, 1117 (1998). 
[5] Y. Rondelez, G. Tresset, T. Nakashima, Y. Kato- 

Yamada, H. Fujita, S. Takeuchi, and H. Noji, Nature 

433, 773 (2005). 
[6] K. Adachi, K. Oiwa, T. Nishizaka, S. Furuike, H. Noji, 

H. Itoh, M. Yoshida, and K. Kinosita, Jr., Cell 130, 309 

(2007). 

[7] E. Muneyuki, T. Watanabe-Nakayama, T. Suzuki, 
M. Yoshida, T. Nishizaka, and H. Noji, Biophys. J. 92, 
1806 (2007). 

[8] T. Harada and S.-I. Sasa, Phys. Rev. Lett. 95, 130602 

(2005) . 

[9] T. Harada and S.-I. Sasa, Phys. Rev. E. 73, 026131 

(2006) . 

[10] J. M. Deutsch and O. Narayan, Phys. Rev. E 74, 026112 
(2006). 

[11] M. Washizu, Y. Kurahashi, H. Iochi, O. Kurosawa, 
S. i. Aizawa, S. Kudo, Y. Magariyama, and H. Hotani, 
IEEE Trans. Industry Appl. 29, 286 (1991). 

[12] R. M. Berry, L. Turner, and H. C. Berg, Biophys. J. 69, 
280 (1995). 

[13] T. Watanabe-Nakayama, S. Toyabe, S. Kudo, 
S. Sugiyama, M. Yoshida, and E. Muneyuki, Biochem. 
Biophys. Res. Comm. 366, 951 (2008). 

[14] C. Bustamante, J. Liphardt, and F. Ritort, Physics To- 
day 58, 43 (2005). 

[15] V. Blickle, T. Speck, L. Helden, U. Seifert, and 
C. Bechinger, Phys. Rev. Lett. 96, 070603 (2006). 

[16] K. Sekimoto, J. Phys. Soc. Jpn. 66, 1234 (1997). 

[17] K. Sekimoto, Stochastic Energetics (Lecture Notes in 
Physics) (Springer, 2010). 

[18] S. Toyabe, H. R. Jiang, T. Nakamura, Y. Murayama, and 
M. Sano, Phys. Rev. E 75, 011122 (2007). 

[19] S. Toyabe and M. Sano, Phys. Rev. E 77, 041403 (2008). 

[20] R. Kubo, M. Toda, and N. Hashitsume, Statsitical 
Physics II (Springer, Berlin, 1991), 2nd ed. 

[21] D. Mizuno, C. Tardin, C. F. Schmidt, and F. C. Mack- 
intosh, Science 315, 370 (2007). 

[22] T. Speck and U. Seifert, Europhys. Lett. 74, 391 (2006). 

[23] V. Blickle, T. Speck, L. Helden, U. Seifert, and 
C. Bechinger, Phys. Rev. Lett. 98, 210601 (2007). 

[24] See EPAPS Document No. [] for Supplementary In- 
formation. For more information on EPAPS, see 
http://www.aip.org/pubservs/epaps.html. 

[25] H. Wang and G. Oster, Europhys. Lett. 57, 134 (2002). 

[26] J. Rosing and E. C. Slater, Biochim. Biophys. Acta 267, 
275 (1972). 

[27] R. W. Guynn and R. L. Veech, J. Biol. Chem. 248, 6966 
(1973). 

[28] O. Panke and B. Rumberg, Biochim. Biophys. Acta 

1322, 183 (1997). 
[29] K. Krab and J. van Wezel, Biochim. Biophys. Acta 1098, 

172 (1992). 



* Electronic address: emuneyuk@phys.chuo-u.ac.jp 



